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ABSTRACT 

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar 
System. These constraints include 1) the orbits, in particular the small eccentricities, and 2) the 
masses of the terrestrial planets - Mars' relatively small mass in particular has not been adequately 
reproduced in previous simulations; 3) the formation timescales of Earth and Mars, as interpreted 
from Hf/W isotopes; 4) the bulk structure of the asteroid belt, in particular the lack of an imprint 
of planetary embryo-sized objects; and 5) Earth's relatively large water content, assuming that it 
was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 
high-resolution (N=1000-2000) dynamical simulations of late-stage planetary accretion with the goal 
of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter 
and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are 
both consistent with and contrary to the "Nice model." We find that a configuration with Jupiter 
and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the 
correct timescale, but Mars' mass is too large by a factor of 5-10 and embryos are often stranded in 
the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly 
higher initial eccentricities (e = 0.07 — 0.1) produces a small Mars, an embryo-free asteroid belt, and 
a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we 
tested reproduced all the observed constraints. Our simulations leave us with a problem: we can 
reasonably satisfy the observed constraints (except for Earth's water) with a configuration of Jupiter 
and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not 
allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which 
are consistent with the Nice model fail to reproduce Mars' small size. 

Subject headings: Terrestrial Planets — Planetary Formation — Accretion — Origin, Solar System 
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1. INTRODUCTION 

It is commonly accepted that rocky planets form by 
the process of coUisional agglomeration of smaller bod- 
ies (for recent reviews, see Chambers 2004, Nagasawa 
et al. 2007 or Raymond 2008). This process starts from 
micron-sized dust grains in young circumstellar disks, 
and the current paradigm proceeds as follows. Grains 
settle to a thin disk midplane on a ~ lO'* year timescale 
(Weidenschilling 1980), and grow quickly via sticky col- 
lisions until they reach cm- or m- sizes (Dullemond & 
Dominik 2004). The time for m-sized bodies to spiral 
in to the star is very short (~ 100 years) such that 
this size range constitutes a barrier to further growth 
(Weidenschilling 1977a). This barrier may be crossed 
by rapid accretion (Weidenschilling & Cuzzi 1993; Benz 
2000) or by local gravitational instability (Goldreich & 
Ward 1973; Youdin & Shu 2002), which can be trig- 
gered by turbulent concentration (Johansen et al. 2007; 
Cuzzi et al. 2008). Larger bodies (100 m to 100 km in 
size) , which are more weakly coupled to the gaseous disk, 
are called planetesimals. Runaway growth of the largest 
planetesimals may occur while the velocity dispersion is 
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small because of strong gravitational focusing such that 
dM/dt - Af4/3 (Safronov 1969; Greenberg et al. 1978). 
However, viscous stirring by the large bodies increases 
the velocity dispersion of planetesimals, thereby reducing 
the growth rate to a roughly geometrical regime, where 
dM/dt - M2/3 (Ida& Makino 1993). Dynamical friction 
acts on the oligarchs, maintaining small eccentricities 
(Ida & Makino 1992; Kokubo & Ida 1998). The building 
blocks of the terrestrial planets, ^Moon-sized planetary 
embryos, form in 10^ — 10® years with a characteristic 
spacing of 5-10 mutual Hill radii (Wetherill & Stewart 
1993; Weidenschilhng et al. 1997; Kokubo & Ida 2000, 
2002). Giant collisions between planetary embryos begin 
to occur when the local density of planetesimals and em- 
bryos is comparable (Wetherill 1985; Kenyon & Bromley 
2006). During late-stage accretion, embryo-planetesimal 
and embryo-embryo impacts are common and the feeding 
zones of terrestrial planets can span several AU in width 
(Raymond et al. 2006a). Late-stage accretion lasts for 
~ 10® years and sets the final bulk architecture of the 
system as well as the composition of the terrestrial plan- 
ets (e.g., Wetherill 1996). 

Past simulations of late-stage accretion have succeeded 
in reproducing several aspects of the Solar System's ter- 
restrial planets. Using only 20-165 particles, Agnor et 
al. (1999) and Chambers (2001) roughly reproduced the 
approximate masses and semimajor axes of Mercury, 
Venus, Earth and Mars. Thommes et al. (2008) also re- 
produced the rough mass distribution of the inner Solar 
System by invoking sweeping secular resonances during 
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the depletion of the Solar Nebula. By taking dynam- 
ical friction from remnant planetesimals into account, 
O'Brien et al. (2006) and Morishima et al. (2008) repro- 
duced the very low eccentricities of the terrestrial plan- 
ets. Several groups have succeeded in delivering water to 
Earth from hydrated asteroidal material, following the 
model of Morbidelli et al. (2000; see also Raymond et 
al. 2004, 2006a, 2007; O'Brien et al. 2006). 

Despite these achievements, no previous study has ade- 
quately reproduced all aspects of the inner Solar System. 
Indeed, as pointed out by Wetherill (1991), Mars' small 
size remains the most difficult constraint to reproduce 
(also discussed in Chambers 2001). Agnor et al. (1999), 
Chambers (2001) and Morishima et al. (2008) succeeded 
in reproducing Mars' small size only because their simu- 
lations started from an annulus of material with a fixed 
width (see also Kominami & Ida 2002). In most cases 
this annulus extended from 0.5-1.5 AU, such that a small 
planet could form at the outer edge of the initial disk be- 
cause of spreading. However, no such edge is thought to 
have existed in the Solar Nebula, so that the assumption 
that the planetesimal and embryo population extended 
only to 1.5 AU is not justified. Chambers (2001) man- 
aged to place Mercury within a planetary mass distri- 
bution but only by adopting an ad-hoc inner disk pro- 
file. Thommes et al. (2008) formed a small Mars but 
the orbits they assumed for Jupiter and Saturn are in- 
consistent with any significant late, planetesimal-driven 
migration of the giant planets (discussed at length in §6.2 
below). In fact, the scenario of Thommes et al. (2008) is 
incompatible with the two currently viable theories for 
the late heavy bombardments because these require ei- 
ther a more compact configuration of Jupiter and Saturn 
(Gomes et al. 2005) or the formation of a small, sub- 
Mars-sized planet at ~ 2 AU (Chambers 2007). 

Terrestrial accretion lasts for ~ 10* years, far longer 
than the few Myr lifetimes of the gaseous component 
of protoplanetary disks (Haisch et al. 2001; Briceno et 
al. 2001; Pascucci et al. 2006). Thus, gas giant planets 
must be fully-formed during late-stage accretion and can 
therefore strongly affect terrestrial bodies, especially if 
the giant planets' orbits are eccentric (Wetherill 1996; 
Chambers & Cassen 2002; Levison & Agnor 2003; Ray- 
mond et al. 2004; O'Brien et al. 2006). Given that suId- 
stantial orbital migration of the Solar System's giant 
planets has been proposed to explain the structure of 
the Kuiper Beh (Fernandez & Ip 1984; Malhotra 1995) 
and the origin of the late heavy bombardment (Strom 
et al. 2005; Gomes et al. 2005), the orbits of Jupiter and 
Saturn at early times are unclear. Indeed, a range of 
Jupiter-Saturn configurations could yield the current So- 
lar System. Thus, if any particular configuration were 
especially adept at reproducing the terrestrial planets, it 
would provide strong circumstantial evidence in favor of 
that configuration. 

In this paper we attempt to reproduce the inner Solar 
System with a suite of high-resolution (N = 1000-2000) 
dynamical simulations of late-stage accretion. We only 
vary one parameter of consequence: the configuration of 
Jupiter and Saturn at early times. We quantify five rele- 
vant constraints that we use to test our models in section 
2. In section 3, we outline our choices of initial condi- 
tions and numerical methods. In section 4 we explore the 
case of two contrasting simulations that each reproduce 



certain constraints. We present results and analysis of 
all simulations in section 5. We discuss these results and 
present our conclusions in section 6. 

2. INNER SOLAR SYSTEM CONSTRAINTS 

We consider five broad attributes which we attempt to 
reproduce statistically with accretion simulations. Other 
observations and measurements exist for inner Solar Sys- 
tem bodies, but we are limiting ourselves to relatively 
broad and well-understood characteristics. These con- 
straints are described below in order from strongest to 
weakest. Weaker constraints rely on models or data that 
are subject to interpretation, while strong constraints are 
directly observed. We use several quantities to compare 
our simulations with the Solar System's terrestrial plan- 
ets. These include statistical measures that were intro- 
duced by Chambers (2001). 

1. The masses and the mass distribution of the terres- 
trial planets. As mentioned above, the mass distri- 
bution of the terrestrial planets, and in particular 
the small masses of Mercury and Mars, have not 
been adequately reproduced in the context of the 
entire Solar System and its history. In this paper 
we do not attempt to reproduce Mercury because 
its small size and large iron content may be the re- 
sult of a mantle-stripping impact (Benz et al. 1988) 
or interesting composition-sorting gaseous effects 
(Weidenschilling 1978). However, for the case of 
Mars, with its more distant orbit, these effects are 
less likely to be a factor, and it should be repro- 
ducible in the context of our simulations. In ad- 
dition, the distribution of mass in the inner Solar 
System is interesting because the majority is con- 
centrated between the orbits of Venus and Earth. 
We therefore use two statistical measures for this 
constraint: 

• The number of planets formed Np. We take 
Np to represent objects that contain at least 
one planetary embryo, that have semimajor 
axes a < 2 AU, and that are on long-term 
stable orbits. It is only for these planets that 
we apply our other measures. 

• A radial mass concentration statistic RMC 
(called in Chambers 2001): 

RMC = max ( ^ J^^j , , (1) 

where Mj and aj are the masses and semi- 
major axes of each planet. The function in 
brackets is calculated for a throughout the 
terrestrial planet zone, and Sc is the maxi- 
mum of that function. This quantity repre- 
sents the degree to which mass is concentrated 
in a small radial annulus: Sc remains small for 
a system of many equal-mass planets but Sc 
is large for systems with few planets and with 
most of the mass in one or two planets. For a 
one planet system, the RMC value is infinite. 
The RMC of the Solar System's terrestrial 
planets is 89.9 (see Table 2). 
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The orbits of the terrestrial planets. The terres- 
trial planets maintain very small orbital eccentric- 
ities and inclinations over long timescales. Earth 
and Venus' time-averaged eccentricities are only 
about 0.03 (e.g., Quinn et al. 1991). Recent simu- 
lations with N > 1000 particles have succeeded in 
reproducing these small eccentricities for the first 
time (O'Brien et al. 2006). We quantify the orbital 
excitation of the terrestrial planets using the nor- 
malized angular momentum deficit AMD (Laskar 
1997). This measures the difference in angular mo- 
mentum of a set of orbits from coplanar, circular 
orbits: 



AMD 



1 — cos{ij)^ 



(2) 



where Oj, e^ , i 



J, and rrij refer to planet j's semi- 



major axis, eccentricity, inclination with respect to 
a fiducial plane, and mass. The AMD of the Solar 
System's terrestrial planets is 0.0018 (see Table 2). 

3. The formation timescales of Earth and Mars. Re- 
cent interpretation of Hf/W measurements sug- 
gest that the last core-formation event on Earth 
occurred at roughly 50-150 Myr (Touboul et 
al. 2007)^5 This event is thought to be the Moon- 
forming impact (Benz et al. 1986; Canup & As- 
phaug 2001). Mars' formation time from Hf/W 
isotopes appears to be significantly shorter, about 
1-10 Myr (Nimmo & Kleine 2007). 

4. The large-scale structure of the asteroid belt. The 
asteroid belt shows a clear division between inner, 
S-types and more distant C-types (e.g., Gradie & 
Tedesco 1982). In addition, there are no large gaps 
in the main belt except those caused by specific 
mean motion or secular resonances. If a planetary 
embryo above a critical mass were stranded in the 
asteroid belt for a long period of time, it would 
disrupt both of these observed characteristics by 
planetesimal scattering (O'Brien et al. , in prepara- 
tion). This constraint puts an upper limit of a few 
lunar masses (~ 0.05 M^) on the mass of an object 
that can survive in the asteroid belt after terres- 
trial planet formation. If an embryo did end up in 
the main belt, it could have been subsequently re- 
moved during the late heavy bombardment (Gomes 
et al. 2005; Strom et al. 2005), but the embryo's dy- 
namical imprint on the asteroid belt would have re- 
mained.^ We note that the asteroid belt is thought 
have been depleted by a factor of ~ 10^ in mass 
over the lifetime of the Solar System. This deple- 
tion is best explained by scattering of planetesi- 
mals by planetary embryos in the primordial belt 

^ Touboul et al. 's (2007) core-formation age is roughly a fac- 
tor of two longer than previous estimates (Kleine et al.2002; Yin 
et al. 2002). It is important to note that Hf/W measurements of 
Earth samples are somewhat uncertain given the unknown amount 
of core/mantle equilibration during giant impacts (Halliday 2004; 
Nimmo & Agnor 2006). However, the samples from Touboul et 
al. (2007) are lunar in origin and therefore circumvent the issue of 
equilibration. 

® It is important to note that the late heavy bombardment 
was a purely dynamical event, as shown by the difference between 
crater size distributions on surfaces older vs. younger than 3.8 Gyr 
(Strom et al. 2005). 
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Fig. 1. — The effect of a Mars-sized planetary embryo on the 
structure of the asteroid belt. Shown are the surviving (massless) 
asteroidal bodies, whose orbits were integrated for 100 Myr under 
the influence of Jupiter and Saturn (not shown). Mars and a Mars- 
mass planetary embryo stranded in the asteroid belt at 2.49 AU. 
Asteroids are color-coded according to their starting semimajor 
axes: grey (2-2.5 AU), light grey (2.5-3 AU), and black (3-3.5 AU). 

(WetheriU 1992; Chambers & WetheriU 2001; Pe- 
tit et al.2001; O'Brien et al.2007), although other 
models do exist (e.g., Lecar & Franklin 1997; Naga- 
sawa et al. 2000). Scattering among embryos often 
places one body in an unstable mean motion reso- 
nances with Jupiter, leading to their rapid removal 
from the belt. This scattering also leads to some 
radial mixing, consistent with the observation that 
the different asteroid taxonomic types are not con- 
fined to narrow zones, but are spread somewhat 
in overlapping but still distinct regions (Gradie & 
Tedesco 1982). Embryos as small as the Moon are 
able to provide the necessary excitation (Chambers 
& WetheriU 2001). Most of the embryos are re- 
moved on a timescale of ~ 10 Myr. However, if 
one or more stray embryos with too large of a mass 
remain in the belt for much longer than this, they 
will lead to excessive radial mixing, inconsistent 
with the observed distribution of different asteroid 
taxonomic types. Figure [1] shows the effect of a 
Mars-mass embryo trapped at 2.5 AU on 100 Myr 
of evolution of 1000 asteroids in the main belt (2- 
3.5 AU), which are assumed to be massless. Two 
features from Fig. [T] are inconsistent with the ob- 
served main belt: the excess radial mixing and the 
gap created in the vicinity of the embryo. More 
massive or more eccentric asteroidal embryos can 
be significantly more disruptive than the case from 
Fig. [U especially if their eccentricity is strongly 
forced by secular perturbations from Jupiter and 
Saturn (O'Brien et al. , in preparation). In addi- 
tion, the simulation from Fig. [1] was only run for 
100 Myr, roughly 500 Myr shorter than the relevant 
timescale, i.e., the time between the completion of 
terrestrial accretion (~ 100 Myr) and the time of 
the late heavy bombardment (600-700 Myr). Thus, 
the constraint we place on our accretion simula- 
tions is that no embryos larger than 0.05 M® can 
survive in the main belt past the end of terrestrial 
planet growth, or in our case 2 x 10^ years. 

5. Earth's water content. One prominent model sug- 
gests that primitive asteroidal material was the 
source of the bulk of Earth's water (Morbidelli 
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et al. 2000; sec also Raymond ct al. 2007). This 
model explains why the D/H ratio of Earth's water 
matches that of carbonaceous chondrites (Robert & 
Epstein 1982; Kerridge 1985), and links Earth's wa- 
ter to the depletion of the primitive asteroid belt. 
Note that other models exist which propose that 
Earth's water came from comets (Delsemme 1994; 
Owen & Bar-Nun 1995), from oxidation of a prim- 
itive, H-rich atmosphere (Ikoma & Genda 2006), 
from adsorption of water onto small grains at 1 AU 
(Muralidharan et al. 2008), or from other sources - 
see Morbidelli et al. (2000) for a discussion of some 
of these models. However, it is our opinion that the 
asteroidal water model of Morbidelli et al. (2000) is 
the most likely source of Earth's water. In fact, wa- 
ter vapor from sublimation of in-spiraling icy bod- 
ies has been detected interior to 1 AU in the pro- 
toplanetary disk around the young star MWC480 
(Eisner 2007); this may be an observation of aster- 
oidal (or in this case potentially cometary) water 
delivery in action. 

3. METHODS 

Our simulations are designed to start at the begin- 
ning of late-stage accretion, after Jupiter and Saturn are 
fully-formed and the nebular gas has dissipated. This is 
probably 1-3 Myr after "time zero." and we base our ini- 
tial conditions on models of the formation of planetary 
embryos (e.g., Kokubo & Ida 2000). We start with a disk 
of planetary embryos and planetesimals, plus Jupiter and 
Saturn. Our simulations are comparable to the highest- 
resolution cases in the literature, containing 85-90 plan- 
etary embryos and 1000-2000 planetesimals.'' 

3.1. Configuration of Jupiter and Saturn 

The resonant structure of the Kuiper Belt appears to 
require a significant outward migration of Neptune (Fer- 
nandez & Ip 1984; Malhotra 1995; Gomes 2003; Levison 
& Morbidelli 2003). This outward migration occurred be- 
cause of the back-reaction from planetesimal scattering, 
which causes the orbits of Saturn, Uranus and Neptune 
to expand and the orbit of Jupiter to contract (Fernandez 
& Ip 1984). In addition, the "Nice model" of giant planet 
evolution, which explains several observed characteristics 
of the Solar System, proposes that Jupiter and Saturn 
formed interior to their mutual 2:1 mean motion reso- 
nance, perhaps in fact in the 3:2 resonance and migrated 
apart (Tsiganis et al. 2005; Gomes et al. 2005; Morbidelli 
et al. 2005, 2007). Thus, Jupiter and Saturn may very 
well have been in a more compact configuration at early 
times. 

We tested a range of configurations for Jupiter and Sat- 
urn, although we did not perform an exhaustive search 
given the large computational expense of each simula- 
tion. However, to account for stochastic variations in 
outcome we performed 4 simulations for each giant planet 
configuration. The configurations we tested were: 

• CJS ("Circular Jupiter and Saturn"). These are 
the initial conditions for the Nice model, as in 

The highest-resolution pubUshed latc-stagc accretion simula- 
tions to date had N=2000-3000 (Raymond et al. 2006a; Morishima 
et al. 2008). 



Tsiganis et al. (2005) and also used in O'Brien et 
al. (2006). Jupiter and Saturn were placed on circu- 
lar orbits with semimajor axes of 5.45 and 8.18 AU 
and a mutual inclination of 0.5 degrees. We note 
that even though Jupiter and Saturn begin with 
zero eccentricities, they induce small, non-zero ec- 
centricities in each others' orbits. 

• CJSECC. Jupiter and Saturn were placed at their 
CJS semimajor axes of 5.45 and 8.18 AU with ej = 
0.02 and eg = 0.03 and a mutual inclination of 0.5 
degrees. 

• E JS ( "Eccentric Jupiter and Saturn" ) . Jupiter and 
Saturn were placed on approximately their current 
orbits: aj = 5.25 AU, ej = 0.05, as = 9.54 AU, 
and es = 0.06, with a mutual inclination of 1.5 

degrees. 

• EEJS ("Extra Eccentric Jupiter and Saturn"). 
Jupiter and Saturn were placed at their current 
semimajor axes but with higher orbital eccentrici- 
ties: aj = 5.25 AU, as = 9.54 AU, and ej = es = 
0.1, with a mutual inclination of 1.5 degrees. These 
cases proved to be interesting, so we ran 8 cases in 
addition to the original four. The next four cases 
(referred to as EEJS 5-8) had the same configura- 
tion of Jupiter and Saturn but 2000 planetesimals 
rather than 1000. The final four cases (EEJS 9-12) 
also had 2000 planetesimals but had ej = 0.07 and 
es = 0.08. 

• JSRES ("Jupiter and Saturn in RESonance"). 

Jupiter and Saturn were placed in their mutual 
3:2 mean motion resonance, following directly from 
simulations of their evolution in the gaseous Solar 
Nebula (Morbidelh et al.2007): aj = 5.43 AU, as 
= 7.30 AU, ej = 0.005, and es = 0.01, with a mu- 
tual inclination of 0.2 degrees. 

• JSRESECC ("Jupiter and Saturn in RESonance 
on ECCentric orbits"). As for JSRES but with 
ej = es = 0.03. 

The EJS and EEJS simulations assume that Jupiter 

and Saturn did not undergo any migration. The EEJS 
simulations are more self-consistent than the EJS simu- 
lations, bec;ause scattering of remnant planetesimals and 
embryos tends to decrease the eccentricities and semi- 
major axes of Jupiter and Saturn (e.g.. Chambers 2001). 
Thus, to end up on their current orbits, Jupiter and 
Saturn would have had to form on more eccentric and 
slightly more distant orbits. The CJS, JSRES and JS- 
RESECC simulations all follow from the Nice model and 
assume that Jupiter and Saturn's orbits changed signifi- 
cantly after their formation, with Saturn migrating out- 
ward and Jupiter inward (Tsiganis et al. 2005) . If mi- 
gration of the giant planets is really associated with the 
late heavy bombardment (Gomes et al. 2005; Strom et 
al. 2005), then at least most of the migration of Jupiter 
and Saturn must have occurred late, well after the com- 
pletion of the terrestrial planet formation process. 

3.2. Properties of the Protoplanetary Disk 
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Fig. 2. — Sample initial conditions for a disk with E ~ r~^l'^ 
containing 97 planetary embryos and 1000 planetesimals. Embryos 
are shown in gray with their sizes proportional to their mass(^/^) 
(but not to scale on the x axis). 

For all of our simulations, the disk of solids extended 
from 0.5 to 4.5 AU and contained populations of plane- 
tary embryos and planetesimals. For most cases, we as- 
sumed that the disk's surface density in solids S followed 
a simple radial power-law distribution: 



(lAu) 



(3) 



For the minimum-mass solar nebula (MMSN) model. 
El « 6 — 7 g cm~^ and x = 3/2 (Weidenschilling 1977a; 
Hayashi 1981). For most of our simulations we assumed 
X = 3/2 but we also performed some cases with a; = 1 for 
the CJS and EJS giant planet configuration. Cases with 
a; = 1 are labeled by the x value; for example, the EJS15 
simulations have x = 3/2 and the EJSl simulations have 
X = 1 (see Table 2). For each case, we calibrated our 
disks to contain a total of 5 in solids between 0.5 
and 4.5 AU, divided equally between the planetesimal 
and embryo components. 

Figure [2] shows a sample set of initial conditions. We 
assumed that embryos are spaced by A = 3-6 mutual Hill 
radii Rh, where Rh = 0.5{ri+r2) [{Mi+ M2)/3Mq]'^^^ , 
where ai and Mi are the radial distance and mass of 
embryo 1. The embryo mass therefore scales with or- 
bital distance as M - ^3/2 (2-^)^3/2 (Xokubo & Ida 
2002; Raymond et al. 2005). The disks contained 85- 
90 embryos with masses between 0.005 and 0.1 M®. In 
Mars' vicinity the typical embryo mass was roughly 1/6 
to 1/3 of a Mars mass. Planetesimals were laid out as 
Np ^ r^+'^ to follow the annular mass, and had masses 
of 0.0025 Mgg. Embryos and planetesimals were given 
randomly-chosen starting eccentricities of less than 0.02 
and inclinations of less than 0.5°. In a few EEJS cases 
we performed additional simulations with 2000 planetes- 
imals, which followed the same distribution but had cor- 
respondingly smaller masses. 

We assume that there existed a radial compositional 
gradient for rocky bodies in the Solar Nebula. This 
gradient was presumably imprinted on planetesimals by 
the local temperature during their formation (e.g.. Boss 
1998), although heating by short-lived radionuclides such 
as ^ Al may have played a role (Grimm & McSween 
1993). We assume the same water distribution as in Ray- 
mond et al. (2004, 2006a), using data for primitive mete- 
orites from Abe et al. (2000). The "water mass fraction", 
WMF, i.e. the water content by mass, varies with radial 



distance r as: 

WMF = 




r < 2 AU 
2 AU < r < 
r > 2.5 AU 



2.5 AU 



(4) 



This water distribution is imprinted on planetesimals 
and embryos at the start of each simulation. During ac- 
cretion the water content of each body is calculated by 
a simple mass balance of all the accreted bodies. We 
do not take into account water loss during giant impacts 
(Genda & Abe 2005; Canup & Pierazzo 2006) or via hy- 
drodynamic escape (Matsui & Abe 1986; Kasting 1988). 

3.3. Numerical Method 

Each simulation was integrated for at least 200 Myr 
using the hybrid symplectic integrator Mercury (Cham- 
bers 1999). We used a 6-day timestep for all integrations; 
numerical tests show that this is adequate to resolve the 
innermost orbits in our simulations and to avoid any sub- 
stantial error buildup (see Ranch & Holman 1999). Col- 
lisions are treated as inelastic mergers, and we assumed 
physical densities of 3gc'm~^ for all embryos and plan- 
etesimals. Simulations were run on individual machines 
in a distributed computing environment, and required 
2-4 months per simulation. The Sun's radius was arti- 
ficially increased to 0.1 AU to avoid numerical error for 
small-perihelion orbits. 

For each Jupiter-Saturn-disk configuration we per- 
formed four different simulations to account for the 
stochastic nature of accretion (e.g.. Chambers & Wether- 
ill 1998). These four cases varied in terms of the random 
number used to initialize our disk code, resulting in dif- 
ferences in the detailed initial distributions of embryos 
and planetesimals. 

Embryo particles interacted gravitationally with all 
other bodies but planetesimal particles did not interact 
with each other. This approximation was made to reduce 
the run time needed per simulation which is already con- 
siderable (see Raymond et al. 2006a for a discussion of 
this issue). The run time r scales with the number of em- 
bryos Ne and the number of planetesimals, Np, roughly 
as T ~ N^ + 2Nf,Np. The non-interaction of planetes- 
imals eliminates an additional N"^ term. Note that r 
refers to the computing time needed for a given timestep. 
The total runtime is r integrated over all timesteps for 
all surviving particles. Thus, a key element in the ac- 
tual runtime of a simulation is the mean particle lifetime. 
Configurations with strong external perturbations (e.g., 
eccentric giant planets) tend to run faster because the 
mean particle lifetime is usually shorter than for config- 
urations with weak external perturbations. 

4. TWO CONTRASTING EXAMPLES 

We illustrate the variations between different cases 
using two simulations with different configurations of 
Jupiter and Saturn: one case from the JSRES batch and 
one from EEJS (simulations JSRES-4 and EEJS-3 in Ta- 
ble 2). Each simulation matched some of our constraints 
but neither matched all of them. Figures [3] and 2] show 
snapshots in the evolution of the two simulations. Prop- 
erties of the planets that formed in each case are listed in 
Table 1. We note that these are individual simulations, 
and that there exists substantial variability in outcome 
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Fig. 3. — Snapshots in time from a simulation with Jupiter and Saturn in 3:2 moan motion resonance (JSRES). The size of each body 
is proportional to its mass'^/'^) (but is not to scale on the x axis). The color of each body corresponds to its water content by mass, from 
red (dry) to blue (5% water). Jupiter is shown as the large black dot; Saturn is not shown. 



TABLE 1 

Planets that formed in the JSRES and EEJS example simulations 



Planet 


a (AU) 


el 


i (deg) 


Mass (M®) 


WMF^ 


Last giant impact (Myr) 


JSRES-a 


0.59 


0.08 


1.7 


0.95 


2.77 X 10-3 


113.5 


JSRES-b 


1.03 


0.03 


2.8 


0.54 


2.87 X 10-2 


160.0 


JSRES-c 


1.42 


0.03 


2.5 


0.85 


5.48 X 10-3 


124.1 


JSRES-d 


1.81 


0.02 


4.7 


0.36 


1.42 X 10-3 


42.9 


EEJS-a 


0.61 


0.08 


3.3 


0.90 


1 X 10-^ 


82.2 


EEJS-b 


1.02 


0.05 


3.2 


0.70 


7.14 X 10-5 


35.6 


EEJS-c 


1.63 


0.16 


9.0 


0.06 


3.08 X 10-2 


0.168 



Orbital values (a, e, i) are averaged over the last 1 Myr of each simulation. 



between simulations even for the same giant planet con- 
figuration. We discuss the outcomes of all simulations in 
section 5. 

In the JSRES simulation (Fig. [3|), eccentricities are 
excited in the inner disk by interactions between em- 
bryos and planetesimals. In the outer disk, eccentricities 
are excited by specific mean motion resonances (MMRs) 
with Jupiter and Saturn: the 3:1, 2:1 and 3:2 MMRs are 
clearly visible. Eccentric embryos perturb nearby bodies 
and act to spread out the resonant excitation on a Myr 
timescale. A stage of chaotic growth lasts for ^100 Myr. 
During this time, there is substantial mixing of objects 
between radial zones, the inner system is cleared of small 
bodies, and four water-rich planets are formed inside 2 
AU with masses between 0.36 and 0.95 (see Table 
1). 



In the EEJS simulation (Fig. S]), the inner and outer 
portions of the disk are quickly divided by a strong secu- 
lar resonance near 2 AU (fg). The evolution of the inner 
disk proceeds in similar fashion to the JSRES simulation, 
although eccentricities are higher because of excitation 
by another secular resonance at 0.7 AU (1^5). The as- 
teroid belt region was cleared more quickly than for the 
JSRES case due to stronger secular and resonant per- 
turbations. The stage of chaotic growth also lasts about 
10® years but with less mixing between radial zones. At 
the end of the simulation, three mainly dry planets have 
formed within 2 AU. The outermost planet lies at 1.63 
AU and is a good Mars analog. 

Figure [5] (top panels) shows the mass of the planets 
over the 200 Myr span of the simulation for the two sim- 
ulations. Planetary growth is a combination of relatively 
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smooth accumulation from a large number of planetesi- 
mals and punctuated accretion from a small number of 
giant impacts with other embryos. In general, embryo- 
embryo collisions increase in magnitude in time simply 
because all embryos are growing. This is particularly 
clear for the case of the innermost planet (0.59 AU) in 
the JSRES simulation which was hit by a 0.41 M® em- 
bryo at 94.8 Myr while the planet was only 0.48 M®. 

The timescale for the last giant impact on the JSRES 
planets was 43-160 Myr, and 0.17-82 Myr for the EEJS 
planets. The Earth analog (i.e., the planet closest to 1 
AU) in each simulation fell slightly out of the 50-150 Myr 
window for the last giant impact on Earth (Touboul et 
al.2007), but on different sides. The JSRES Earth ana- 
log's last giant impact was slightly too late (160 Myr) 
while the EEJS Earth analog's was too early (35.6 Myr). 
The Mars analog in the JSRES simulation (at 1.42 AU) 
has a mass that is roughly eight times too large and a 
formation timescale that is far too long (124 Myr as com- 
pared with the Hf/W isotopic age of 1-10 Myr; Nimmo 
& Kleine 2007). In contrast, the EEJS simulation pro- 
duced an excellent Mars analog that is actually somewhat 
smaller than Mars (0.06 M® vs. 0.11 M®) and whose 
only giant impact occurred 168,000 years into the simu- 
lation. Given that the "time zero" for our simulations is 
probably 1-3 Myr after the formation of the Solar Neb- 
ula, this is consistent with isotopic measurements. It is 
interesting to note that the last giant impact on the in- 
nermost planet in each simulation occurred quite late, at 
~ 10^ years (see Table 1). The reason for the late im- 
pact was different for the two simulations. For the JSRES 
simulation the last giant impactor originated in the as- 



teroid belt, where the timescale for close encounters and 
scattering is longer than the inner system. For the EEJS 
simulation, the last giant impactor originated at 1.2 AU 
but had its inclination increased by a short time spent 
in the vicinity of the secular resonance, thereby pro- 
longing its dynamical lifetime in the inner system. These 
late giant impacts on close-in planets contrast with the 
nominal view of accretion occurring fastest in the inner 
regions of the disk, especially given the much shorter ac- 
cretion timescales for the Earth and Mars analogs in the 
EEJS simulation. 

The feeding zones of the planets from the JSRES and 
EEJS simulations are shown in the bottom panels of 
Fig. m Feeding zones were calculated as the fraction 
of material incorporated into each planet that originated 
in each 0.45 AU-wide radial bin. The feeding zones of all 
planets overlap in each simulation, although the width 
of individual feeding zones vary.* In the JSRES simu- 
lation each of the four planets accreted material from a 
radial width of more than 3 AU. Given that the source 
of water lies beyond 2-2.5 AU, this explains the large 
water abundance in the JSRES planets. The Earth ana- 
log's feeding zone is exceedingly wide and is unusual in 
that its accretion seed actually started the simulation in 
the outer asteroid belt, at 4.3 AU.^ The Earth analog's 
water content was therefore very large, roughly 30 times 
the Earth's current water content without accounting for 

* Terrestrial feeding zones are not static, but actually widen and 
move outward in time (Raymond et al. 2006a). 

^ A planet's accretion seed is simply the object that was the 
larger in each of its collisions. The planet retains the name of this 
object. 




Fig. 5. — Growth and feeding zones for the planets that formed in our example JSRES and EEJS simulations. The top panels show the 
growth of each surviving planet interior to 2 AU for the JSRES (left) and EEJS (simulations). The 50-150 Myr isotopic constraints on the 
timing of the Moon-forming impact on Earth are shaded. Each colored curve corresponds to an individual planet, as labeled. The bottom 
panels show the origin of the material incorporated into the planets. See Table 1 and the text for details. 



any water loss (the Earth's WMF is ^ 10^^; Lecuyer et 
al. 1998). In contrast, the three planets from the EEJS 
shTiulations each had feeding zones of less than 1.7 AU 
in width. Very little material from exterior to 2 AU was 
incorporated into the EEJS planets, with the notable ex- 
ception of one embryo that originated at 2.64 AU and was 
the accretion seed of the Mars analog. Thus, the Earth 
and Venus analogs are very dry, but the Mars analog is 
very water-rich. 

The JSRES and EEJS simulations each reproduced 
some of our constraints but neither reproduced them all. 
The JSRES simulation formed a terrestrial planet sys- 
tem with eccentricities and inclinations almost as low as 
the Solar System's terrestrial planets', with an AMD 
of 0.0023 (as compared with 0.0018 for Mercury, Venus, 
Earth and Mars - hereafter MVEM). The JSRES plan- 
ets also contain abundant water that was delivered from 
the primordial asteroid belt. The formation timescale 
of the Earth analog is roughly consistent with isotopic 
constraints for Earth. However, the JSRES Mars analog 
bears little resemblance to the real planet in terms of its 
mass and formation timescale. In addition, three extra 
large bodies exist at the end of the simulation: a 0.36 
M0 planet at 1.8 AU and two embryos in the asteroid 
belt totaling 0.11 Mq. These remnant bodies, in par- 
ticular the embryos in the asteroid belt, are inconsistent 
with the observed inner Solar System. 

The EEJS planets are novel among accretion simu- 
lations of this kind because they contain a reasonable 
Mars analog in terms of its mass, orbit and formation 
timescale. In addition, the approximate masses and spac- 



ing of the EEJS planets are close to those of Venus, Earth 
and Mars. No embryos are stranded in the asteroid belt 
although a dozen planetesimals remain in the belt. How- 
ever, the EEJS Earth analog's formation timescale is too 
short by ^ 20%. More importantly, the AMD for the 
system is 0.0086, roughly 5 times higher than the MVEM 
value. Finally, the Earth analog is almost completely de- 
void of asteroidal water and thus requires an alternate 
source. 

The values of the radial mass concentration statistic 
RMC of both simulations are far lower than for the inner 
Solar System. The RMC values are 28.5 for the JSRES 
simulation and 44.2 for the EEJS case, as compared with 
89.9 for MVEM. 

Figure [6] shows the details of each planetesimal and 
embryo collision that occurred on the surviving planets 
in the two simulations. The impact angle theta is defined 
to be zero for a head-on collision and 90° for a grazing 
collision. The impact velocity is given in terms of the 
two-body escape speed v^sc ■ 

he {M1 + M2) 

where G is the gravitational constant. Mi and M2 are the 
colliding bodies' masses, and i?i and R2 are the bodies' 
physical radii. In the absence of 3-body effects, which are 
relevant in at most a few percent of collisions, collisions 
can only occur at v/vgsc > 1- 

Planetesimal-embryo impacts (small symbols in Fig. [S]) 
tend to occur at higher velocities than embryo-embryo 
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Fig. 6. — Impact angles and velocities for all collisions that 
occurred in our example EEJS and JSRES simulations. Large grey 
symbols refer to embryo-embryo impacts and small black symbols 
to embryo-planetesimal impacts. Each symbol refers to the impacts 
that occurred on a specific planet - for each simulation, planets are 
ordered by their proximity to the star: planet a is closest, followed 
by planet b, etc (the EEJS simulation only formed three planets 
so there is no planet d). Note that theta=0 and 90° corresponds 
to head-on and graxing impacts, respectively. Impact velocities 
are normalized by the two-body escape speed Vesc (see Eqn 5). 
The dashed line is the approximate boundary between accretionary 
(below the line) and erosive collisions (Agnor & Asphaug 2004). 

impacts simply because their eccentricities are higher on 
average due to viscous stirring. However, we expect vir- 
tually all planetesimal-embryo collisions to result in net 
growth. Note that our numerical scheme does not allow 
for planetesimal-planetesimal collisions (see discussion in 
§6). 

High-speed or off-center embryo-embryo collisions 
(large symbols in Fig. [S]) can result in either partial accre- 
tion or even erosion. Agnor & Asphaug (2004) showed 
that accretionary collisions only occur at v/vesc ^ 1-5 
and preferentially for small impact angles. The majority 
of giant impacts (filled circles in Fig. [6]) occur at low 
speeds and should therefore be accretionary. Indeed, 
the majority of impacts lie in the accretionary regime 
as defined by Agnor & Asphaug (2004): 60%, 71%, 71%, 
and 75% for the four planets from the JSRES simulation 
(listed from closest- to farthest from the Sun), and 83%, 
100% and 0% for the three EEJS planets. This is a 
larger fraction than the 55% found by Agnor et al. (1999) 
and the ~ half inferred by Agnor & Asphaug (2004). 
We assume that dynamical friction from small bodies re- 
duced the mean impact speed and increased the fraction 
of accretionary impacts. We note, however, that the def- 
inition of an accretionary impact only requires that the 
collision produce an object larger than either of the two 
impactors, not that the object's mass equal the sum of 
the colliding masses. In particular, each of our example 
simulations has a cluster of impacts at low velocity but 
large angle (see discussion in section 5). A large number 
of fragments is probably produced in these off-center col- 
lisions (Asphaug et al. 2006). We do not have the ability 
to track the effect of these fragments, which could be 
important (see discussion in Section 6). On the other 
hand, an erosive, "hit and run" collision usually results 
in an extra body that looks very similar to the original 

Note that for the outermost EEJS planet (the Mars analog), 
only one giant collision occurred, when the planet was just 0.036 

and was hit by a 0.015 Mgj embryo. The collision was nearly 
head-on (sintheta = 0.14) but high-speed {v/vesc = 5.2). 
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Fig. 7. — Orbital eccentricities of massless test particles in the 
inner Solar System after 1 Myr for the JSRES (left panel) and 
EEJS (right panel) configurations of Jupiter and Saturn. Note the 
difference in the y axis scale between the two panels. 

impactor and can easily be accreted in a later collision 
involving that extra body (Asphaug et al. 2006). 

The different outcomes in the EEJS and JSRES exam- 
ple simulations can be attributed to differences in eccen- 
tricity and inclination excitation by specific resonances 
with Jupiter and Saturn, as well as by secular pertur- 
bations from Jupiter and Saturn. Figure [7] shows the 
eccentricities of test particles on initially circular orbits 
after 1 Myr of evolution in each giant planet system, with 
no embryos present. For the JSRES case, the amount 
of eccentricity excitation is small. The main sources of 
excitation are the secular resonance at 1.3 AU and 
the 2:1 MMR with Jupiter at 3.4 AU. The 3:1 and 3:2 
MMRs with Jupiter are faintly visible at 2.6 AU and 4.1 
AU. The small amount of external forcing means that 
the self-scattering of embryos and planetesimals is the 
dominant source of eccentricity in the JSRES simula- 
tions. Given that the disk is continuous and contains a 
significant amount of mass in the Mars region, no dynam- 
ical mechanism exists to remove that mass. In addition, 
the weak influence of the giant planets allows for effi- 
cient delivery of water-rich material via a large number of 
relatively weak embryo-embryo and embryo-planetesimal 
scattering events (Raymond et al. 2007). 

In contrast, perturbations from Jupiter and Saturn 
play a dominant role in the EEJS configuration. Strong 
secular resonances are visible at 0.7 AU (v^) and 2.2 AU 
(z/g). In addition, secular excitation is strong enough to 
impart a typical free eccentricity of 0.1-0.2 throughout 
the inner Solar System. The i>e secular resonance is di- 
rectly responsible for Mars' small size, as it efficiently re- 
moves mass from the 1.5-2.5 AU region, mainly by driv- 
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ing eccentricities of bodies to 1 and inducing collisions 
with the Sun. However, the vq acts as a barrier between 
the terrestrial planets and the asteroid belt, such that 
water delivery is severely reduced. The strong eccentric- 
ity forcing throughout the inner Solar System appears 
to prevent low-AMD terrestrial planets from forming. 
However, the scattering of embryos and planetesimals by 
Jupiter and Saturn throughout accretion reduces the gi- 
ant planets' eccentricities and weakens their secular per- 
turbations in time. 

Despite the differences between the JSRES and EEJS 
simulations, it is important to realize that small secular 
perturbations from the giant planets do not necessarily 
correlate with low- AMD terrestrial planets, especially in 
the case of limited numerical resolution. In fact, O'Brien 
et al. (2006) formed significantly lower-^MI? terrestrial 
planets for the EJS configuration than their CJS sim- 
ulations. The reason for this is that the timescale for 
the removal of asteroidal material was very long for the 
CJS simulations, and close encounters with late-arriving 
material from the asteroid belt tended to increase eccen- 
tricities. In contrast, their EJS simulations cleared out 
the asteroid belt quickly and the secular forcing of eccen- 
tricities was small enough inside ~ 2 AU that the giant 
planets did not act to increase the terrestrial planets' 
AMD. 

5. SIMULATION OUTCOMES AND COMPARISON WITH 
OUR CONSTRAINTS 

The evolution of each simulation proceeded in a qual- 
itatively similar fashion to the example EEJS or JSRES 
simulations. In fact, the two cases illustrated in section 4 
comprise the most extreme variations in our sample. The 
other cases lie between those extremes, typically with a 
moderate amount of excitation from the giant planets in 
the outer disk and relatively little external excitation in 
the inner disk. In this section we discuss the outcomes 
of our simulations in terms of how they compare with 
our Solar System constraints. We explore how the differ- 
ences between cases can be attributed to the giant planet 
configuration and, to a lesser degree, to variations in the 
disk's density profile. 

There was a large range in the characteristics of the 
terrestrial planet systems that formed. The number of 
planets in a given system ranged from 2-6, where we 
define a planet to contain an least one embryo, to be 
interior to 2 AU, and to be on a stable orbit that does 
not cross the orbit of any other planets or embryos. The 
total mass in planets varied by almost a factor of two, 
from 1.4 to 2.7 M®. Table 2 summarizes the outcome of 
each simulation. 

Figure [8] shows the median collision velocities and an- 
gles for the giant (embryo-embryo) collisions that oc- 
curred during the formation of the Earth and Mars 
analogs in each of our simulations."'^^ Despite the ex- 
istence of high- velocity impacts (see Fig. the me- 
dian coUisional values are quite modest and in almost 
all cases the vast majority of giant collisions are accre- 
tionary rather than erosive. 

Earth and Mars analogs are defined to be the most massive 
planets in the region from 0.8-1.25 AU, and 1.25-1.75 AU, respec- 
tively. If no planet exists in that zone, then the Earth analog is 
taken to be the planet that is closest to 1 AU and the Mars analog 
is taken to be the outermost planet inside 2 AU. 
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Fig. 8. — Median angles and velocities for the giant collisions 
that formed the Earth and Mars analogs in each simulation, la- 
beled with different symbols (see legend). The shaded region rep- 
resents the zone where impacts should be accretionary rather than 
erosive (Agnor & Asphaug 2004). Recall that theta=0 and 90° 
corresponds to head-on and graxing impacts, respectively. Impact 
velocities are normalized by the two-body escape speed Vesc (see 
Eqn 5). 

If typical impact speeds or angles on Mars analogs were 
much higher than for Earth analogs, then the fraction 
of erosive collisions on Mars analogs would be higher 
and one could claim that the mass ratio of Mars- to 
Earth-analogs in our simulations was too high. How- 
ever, the distribution of collision velocities for Earth and 
Mars analogs is very similar (Fig. [5]). Therefore, we can 
rule out variations in impact properties as the source of 
large Mars analogs in our simulations. One exception 
are the EJS simulations, many of which have somewhat 
higher impact speeds and angles for Mars analogs than 
for Earth analogs. For the EJS simulations, the mass of 
Mars analogs may therefore be somewhat overestimated. 

Canup (2004) showed that a very particular impact 
configuration was required to form the Moon. Such an 
impact must be low- velocity {v/vesc < 1-1), off-center 
(sin[theta] between 0.67 and 0.76), and have an impactor 
to target mass ratio between 0.11 and 0.15. Canup (2008) 
found that, for prograde rotation of the proto-Earth, 
slightly smaller impactors can form Moon-analogs, with 
a cutoff at roughly 0.1. For retrograde rotation of the 
proto-Earth, larger impactor-to-target mass ratios are al- 
lowed but ratios less than 0.1 are still unable to form the 
Moon. We examined the last three giant impacts suffered 
by the Earth analog in each simulation, and none of the 
impacts fulfilled Canup's (2004) three requirements. In 
fact, none of the last three impacts on an Earth analog 
had an impactor to target mass ratio larger than 0.05 for 
the right collision angle and speed. In addition, only 4% 
of the late giant impacts satisfied Canup's velocity and 
angle criteria. We conclude that the Earth's Moon must 
be a cosmic rarity unless differences between planetary 
systems produce a systematic change in the likelihood of 
Moon-forming impacts. However, we note that the simu- 
lations of Canup (2004, 2008) were specifically designed 
to reproduce the details of the Earth-Moon system, in 
particular its high specific angular momentum and the 
Moon's small core. We suspect that a much larger range 
of late giant collisions would produce satellites, although 
their properties could be much different than the Moon. 

As noted in section 4, our simulations show a group- 
ing of low-velocity grazing collisions. The difference be- 
tween the velocity distributions at small and large angles 
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Comparison between simulations and observed constraints^ 
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JSRES-l 


6 


2.61 


0.0022 


32.4 


7.5 X 10-3 


0.54 


29.2 


28.1 


2 


JSRES-2 


2 


2.31 


0.0119 


61.0 


5.3 X 10-3 


1.27 


176.8 


176.8 


6 


JSRES-3 


4 


2.55 


0.0071 


34.2 


1.2 X 10-3 




115.3 




1 


JSRES-4 


4 


2.70 


0.0023 


28.5 


2.9 X 10-2 


0.85 


160.1 


124.1 


2 


JSRESECC-1 


4 


2.70 


0.0016 


28.7 


1.3 X 10-3 


1.01 


20.2 


81.0 


1 


JSRESECC-2 


4 


2.73 


0.0044 


26.8 


3.5 X 10-4 


0.96 


73.7 


60.6 


1 


.JSRESECC-3 


4 


2.51 


0.0041 


28.1 


3.0 X 10"" 


0.73 


78.5 


99.8 


3 


JSRESECXJ-i 


3 


2.18 


0.0025 


:^<).:^ 


1.1 X 10 


0.9S 


110.3 


177.9 


2 


MVEM3 


4 


1.98 


0.0018 


89.9 


~ 1 X 10-3 


0.11 


50-150 


1-10 






"Table columns arc: the simulation, the number of terrestrial planets inside 2 AU Npi, the total mass in those planets Mtot, the angular 
momentum defieit AMD (see Eqn 2), the radial mass coneentration statistic RMC (see Eqn 1), the water content by mass of the simulation's 
Earth analog WMF(^, the mass of the simulation's Mars analog, the time of the last giant impact on the Earth and Mars analogs T/q 
and TfM, and the number of embryos stranded in the asteroid belt that were more massive than 0.03 Mq. A comparison with the Solar 
System's terrestrial planets (MVEM) is shown at the bottom. 

''A few simulations did not form Mars analogs at all (i.e., no planets between 1.25 and 1.75 AU). 

"^Earth's water content is not well known because the amount of water in the mantle has been estimated to be between 1-10 "oceans", 
where 1 ocean (= 1.5 x lO^^p) is the amount of water on Earth's surface (Lecuyer et al. 1998). Our estimate of 10-3 for Earth's water 
content by mass assumes that 3 oceans are locked in the mantle. 

is only significant for grazing collisions with sin (theta) 
> 0.9, where the collisions that were registered by the 
code do indeed occur at lower speeds. When comparing 
the statistics of head-on (sin [theta] < 0.7) and grazing 
(sin [theta] > 0.9) collisions in all 40 simulations, there 
were no notable differences in terms of collision time, 
distance from the Sun, or the details of the impactor. 
Grazing collisions did, however, have target masses ~ 
15% larger than for head-on collisions. Statistically, one 
would expect collisions between equal-mass objects to 
have a wider distribution in sin(theta) than for collisions 
with one dominant mass, although this also depends on 
the collision speed. The small mass ratio for the graz- 



ing collisions may explain the low collision speeds, sim- 
ply because the two-body escape speed is larger than for 
unequal-mass objects. 

It is possible that the Mercury code (Chambers 1999) 
has difficulty registering high-speed grazing collisions 
because they could travel many Hill radii in a single 
timestep. For example, two bodies traveling with rel- 
ative velocity of 10 kms^^ travel 0.035 AU with respect 
to each other in a single 6 day timestep. The Earth's Hill 
sphere Rr is about 0.01 AU, and approaches within 2>Rh 
are tracked numerically with the Bulirsch-Stoer method 
rather than the symplectic map. Thus, any two objects 
that are flagged as having a close encounter will have 
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their orbits faithfully resolved. However, if the two ob- 
jects were never flagged to approach within 3Rh then an 
encounter could be missed. If that were the case then 
we would expect to miss more collisions at small orbital 
distances because the Hill sphere is smaller and relative 
velocities are larger. Although our statistics are limited, 
we don't see any evidence for this. On the other hand, 
the easiest grazing collisions for Mercury to find should 
be those between massive bodies traveling at low speeds, 
and we can think of no obvious physical reason that high- 
speed grazing collisions should not occur. Thus, although 
we have not found any evidence of the code missing high- 
speed grazing collisions, we can not rule out the possibil- 
ity. We expect that such collisions would likely result in 
a "bounce" rather than a collision (Asphaug et al. 2006), 
and that later lower-speed or head-on collisions, impacts 
that are certainly found by Mercury, would be the ones 
to result in planetary growth. 

The effect of varying the disk's surface density pro- 
file between (i.e., x = 1 from Eqn. 3) and r~^-^ 
(x = 3/2) not insignificant, although we only varied this 
parameter for the CJS and EJS configurations of Jupiter 
and Saturn (EJSl and CJSl had r"!, EJS15 and CJS15 
had r-i-5). The properties of the CJSl, CJS15, EJSl and 
EJS15 simulations are summarized in Table 3. The 
simulations formed slightly fewer planets, contained less 
total mass in planets, had longer formation timescales for 
Earth and higher AMD values for the final systems than 
the r~^/^ simulations. In addition. Earth analogs in the 
EJSl simulations contained far more water than Earths 
in the EJS15 simulations.^^ These trends are consistent 
with the results of Raymond et al. (2005), and appear to 
be due simply to the fact that the simulations con- 
tain far more mass in the asteroid region than the r~^-^ 
simulations. Given that planets form largely from local 
material, the simulations contain less material in the 
inner disk and therefore form less massive planets. In ad- 
dition, the relatively large amount of water-rich material 
in the asteroid belt increases the probability of water de- 
livery, although water delivery is more sensitive to the 
giant planet configuration than to the disk properties: 
despite the large median value, note that 2/4 EJSl simu- 
lations formed Earths with less than 1 part per thousand 
of water (see Table 2). The large amount of asteroidal 
material in the simulations also prolongs the period 
of chaotic bombardment, increasing the mean formation 
time for Earth. 

The higher AMD values for r^^ simulations appears 
to be linked to the mean formation timescale. Indeed, 
Figure [5] shows a weak correlation between the timescale 
for the last giant impact on Earth Tform,® ^-nd the AMD 
of the system for all 40 simulations. For Tform,® < 100 
Myr, the median AMD is 0.003 and for Tform]® > 100 
Myr, the median AMD is 0.010. We attribute this trend 
to the fact that the planetesimal population decays with 
time and our simulations have limited resolution. So for 
late giant impacts or scattering events among the em- 
bryos, there are fewer planetesimals around to re-damp 

The same trend was not seen in the CJSl vs. CJS15 simula- 
tions because the small number of planets that formed in the CJSl 
cases led to the Earth analog being located at ~ 0.8 AU in 3/4 
cases. Given that water delivery decreases with distance from the 
water source (> 2 — 2.5 AU), this decreased the water content of 
Earth analogs in the CJSl simulations. 
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Fig. 9. — The system angular momentum deficit, normalized to 
the MVEM value of 0.0018, as a function of the time of the last 
giant impact for the Earth analog in each of our simulations. The 
symbol for each simulation is the same as in Fig. |8] The region of 
successful outcomes is shaded. 

the planets by dynamical friction if the planets form more 
slowly. This problem could be alleviated with simula- 
tions which continuously regenerate planetesimals from 
the debris of giant impacts (e.g., Levison et al. 2005), be- 
cause in that case the planetesimal population would be 
sustained for as long as the giant impacts occur. 

For the remainder of our analysis we consider the gi- 
ant planet configuration as the only variable. We there- 
fore combine the CJSl and CJS15 simulations into CJS, 
and the EJSl and EJS15 cases into EJS. Given that the 
failings of these simulations are generally the same (see 
Table 2), we do not expect this to skew our results. 

Figure [TO] shows the mass-semimajor axis distribution 
of all 40 of our simulations, grouped into categories 
with similar giant planet configurations. When com- 
pared with the Solar System's terrestrial planets, it is 
immediately evident that all cases with circular or lower- 
eccentricity giant planets fail miserably at reproducing 
Mars' small size. Indeed, for the CJS, CJSECC, JS- 
RES and JSRESECC simulations, planets in Mars' vicin- 
ity are typically 0.5-1 M^. The radial distribution of 
massive planets is much broader in these cases than in 
the Solar System, and Earth-sized planets are commonly 
formed all the way out to 2 AU. In contrast, the EJS and 
especially the EEJS simulations did a much better job of 
reproducing Mars' small size. The EJS simulations have 
smaller Mars analogs than the CJS and JSRES cases but 
in most cases MMars ~ 0.3 M^. In 5/12 EEJS simula- 
tions the Mars analog was between 0.06 and 0.2 M^, 
and in each of those 5 cases the last giant impact oc- 
curred before 10 Myr. The radial mass distributions for 
the EJS and EEJS simulations are peaked, as is the case 
for MVEM. For the EJS simulations the peak is close to 
(or perhaps slightly interior to) 1 AU, but for the EEJS 
simulations the most massive planets tend to lie inte- 
rior to 1 AU and planets at 1 AU are typically half an 
Earth mass. This can be explained as a byproduct of 
the excitation of the planetesimals and embryos by the 
giant planets: planetesimals and embryos on eccentric 
orbits are most likely to collide close to perihelion, such 
that systems with eccentric giant planets tend to have 
the most massive planets closer to their stars than for 
systems with low-mass or low-eccentricity giant planets 
(Levison & Agnor 2003). 

It is also clear from Fig. [TO] that simulations with cir- 
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TABLE 3 

Mean properties of terresrtial planet systems for different disk surface density profiles-' 



Configuration 



Mean Nr, 



Mean Mtot(Mi^) 



Median AMD 



Median VKMFq 



Median T, 



f orm.i 



)(Myr) 



CJSl 
CJS15 

EJSl 
EJS15 



2.75 
3.5 

2.5 
3.0 



2.32 
2.77 

1.55 
2.06 



0.017 
0.003 

0.010 
0.005 



3.2 X 10-3 
5.7 X 10-3 

6.3 X 10-3 
2.6 X 10-* 



140 
104 

169 
66 



"Recall that the CJSl and EJSl sims had disks with r ^ surface density profiles, while CJS15, EJS15, and all our other simulations had 
r~^-^ surface density profiles. 
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Fig. 10. — Mass vs. semimajor axis for a range of simulations with different configurations of Jupiter and Saturn. Each panel shows all 
planets that formed in the relevant simulations (see Table 2) as grey circles, with horizontal lines representing the orbital eccentricity. The 
Solar System's terrestrial planets are shows as the black squares, with 3 Myr averages for their eccentricities in grey (taken from Quinn et 
al. 1991). 



cular or low-eccentricity giant planets (CJS, CJSECC, 
JSRES, JSRESECC) tend to strand massive embryos in 
the asteroid belt. These embryos are typically 0.05-0.2 
and would certainly disrupt the observed asteroid 
distribution. In contrast, the EJS and EEJS simulations 
leave fewer embryos in the asteroid belt, and those that 
are stranded are typically smaller. 

A trend that is less evident from Fig. [10] is that the 
total terrestrial planet mass decreases with the giant 
planet eccentricity in almost all cases. The one excep- 
tion to this rule is for the JSRESECC simulations, which 
have roughly the same total planet mass as the JSRES 
cases; however, at the end of each JSRESECC simulation 
ej < 0.01, so the difference between the two giant planet 
configurations is actually fairly minor. The reason for the 
correlation between increased giant planet eccentricity 
and decreased total mass in terrestrial planets is simply 
that eccentric giant planets perturb terrestrial and aster- 



oidal bodies more strongly and destroy a larger fraction 
of the disk via ejection and collisions with the Sun than 
circrdar giant planets (Chambers & Cassen 2002; Levison 
& Agnor 2003; Raymond et al. 2004; O'Brien et al. 2006). 
In addition, terrestrial planets can't form as close to ec- 
centric giant planets as they can to circular giant planets 
(Raymond 2006); this may explain the reduced number 
of stranded asteroidal embryos for the EJS and EEJS 
simulations. 

Figure [11] shows the radial mass concentration statis- 
tic RMC as a function of the angular momentum deficit 
AMD for the terrestrial planet system that formed in 
each simulation. These statistics are normalized with 
respect to the MVEM values of 0.0018 and 89.9. The 
systems that formed have a wide range in AMD, from 
0.5 to almost 20 times the MVEM value. In contrast, 
systems are clumped in RMC between 0.3 and 0.7 times 
the MVEM value; none has RMC higher than 0.71 - this 
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Fig. 11. — The angular momentum deficit AMD (Eqn. 2) vs. 
radial mass concentration statistic RMC (Eqn. 1) for all of our 
simulations. Both the AMD and RMC are normalized to the 
MVEM values of 0.0018 {AMD) and 89.9 {RMC). The terrestrial 
planet system from each simulation is represented as a single point, 
following the label at the bottom of the plot (e.g., EEJS simulations 
are filled circles, etc). The region of successful outcomes is shaded. 

is similar to the results of Chambers (2001). Most sys- 
tems have AMD values somewhat larger than MVEM, 
although a few cases have AMD smaller than MVEM 
(see Table 2). 

Table 4 lists the median AMD and RMC values for 
each giant planet configuration. The AMD varies sig- 
nificantly for the different configurations, and is smallest 
for the JSRESECC and CJSECC simulations. This is 
consistent with the results of O'Brien et al. (2006), who 
formed lower- AMD systems for cases with moderately 
eccentric giant planets (analogous to our EJS simula- 
tions). The most eccentric planets formed in the CJS 
and EEJS simulations. The CJS systems were eccen- 
tric because the timescale for clearing out of the aster- 
oid belt and late encounters with remnant embryos was 
relatively long such that few planetesimals remained for 
damping after late encounters, as discussed above. The 
EEJS systems were eccentric mainly because of the ex- 
citation caused by the v^, and v^, resonances, as well as 
direct perturbations by the giant planets. However, our 
numerical resolution also played a significant role: the 
four EEJS simulations with 1000 planetesimals (EEJS 
1-4) had a median AMD of 0.015, but the four simula- 
tions with 2000 planetesimals and slightly less eccentric 
giant planets (EEJS 9-12) had a median AMD of 0.0033. 
Note, however, that the EEJS terrestrial planet systems 
with 2000 planetesimals and ej = eg = 0.1 (EEJS 5- 
8) were even slightly more eccentric than the cases with 
1000 planetesimals. This large variation in AMD for the 
EEJS cases is again linked to the presence or absence of 
damping at the time of the last encounters between em- 
bryos. Indeed, the mean formation timescale for Earth 
analogs again scales inversely with the AMD: for EEJS 
1-4, 5-8, and 9-12 the median formation timescales for 
Earth were 109, 129, and 55 Myr. This continues an im- 
portant trend that we see: longer formation timescales 
lead to higher AMD because fewer planetesimals exist 
for dynamical friction at later times. This trend is caused 
in part by our numerical resolution (1000-2000 planetes- 
imals instead of billions) and in part because we do not 
account for impact debris. However, the planetesimal 
population certainly does contribute to decreasing eccen- 
tricities so we believe that this effect is real, although ac- 



FlG. 12. — The water content by mass W MF^ vs. the time of the 
last giant impact for the Earth analog in each of our simulations. 
The sy mb ol for each simulation is the same as in Fig. 1111 Also as 
in Fig. 111! the region of successful outcomes is shaded. 

counting for other factors should weaken the correlation. 

We performed a suite of two-sided Kolmogorov- 
Smirnov and Wilcoxon tests to determine which differ- 
ences in AMD and RMC in our simulations were sta- 
tistically significant. At the 0.05 level (i.e., p < 0.05 
where p is the probability that the two distributions were 
drawn from the same sample), we only found differences 
between the EEJS normal-resolution simulations (EEJS 
1-4 = EEJSnr [with 1000 planetesimals], EEJS 5-12 = 
EEJShr [2000 planetesimals]) and CJSECC, EJS15, and 
JSRESECC, all with p < 0.029 from Wilcoxon tests. For 
the RMC values, the following sets of simulations pro- 
vided significantly different values {p < 0.05): CJS15 
vs. EJSl ip < 0.029), CJS15 vs. EEJShr {p < 0.016), 
EEJShr vs. JSRESECC [p < 0.016), and EJSl vs. JS- 
RESECC {p < 0.029). If we only include variations of 
giant planet configuration and ignore changes in disk 
surface density profile (for the CJS and EJS simula- 
tions) and resolution (for the EEJS simulations), there 
are no statistically significant differences in AMD, but 
for RMC there are differences between the following con- 
figurations: CJS vs. CJSECC (p < 0.048), CJS vs. 
EJS {p < 0.007), CJS vs. EEJS {p < 0.005), EJS vs. 
JSRESECC {p < 0.016), and EEJS vs. JSRESECC 
{p < 0.008). Thus, the normal- resolution EEJS simu- 
lations have significantly higher AMD values than the 
other simulations, but increasing the resolution brings 
them into agreement with the other cases (just as in- 
creasing the resolution for the other cases would likely 
also decrease their AMD values). The CJS (especially 
CJS 15) and JSRESECC simulations represent the sta- 
tistically smallest RMC values of our sample, and the 
highest come from the EJS and EEJS simulations, al- 
though these arc still far below the MVEM value. 

We follow Thommes et al. (2008) and define a success- 
ful outcome as a system for which both the AMD and 
RMC are within a factor of two of the MVEM values. 
This successful area is shaded in Fig. [TT] Only four sim- 
ulations were successful in terms of AMD and RMC , 
and all had relatively eccentric giant planets: EJSl-4, 
EJS15-2, EEJS-U and EEJS-12. 

Figure \n\ shows the water content by mass WMF(^ of 
the Earth analog in each simulation vs. the time of the 
last giant impact on that same planet Tform,®- There 
is a wide range in both of these parameters, and some 
correlation with the giant planet configuration (see Table 
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TABLE 4 

Statistical values of terrestrial planet systems for different giant planet configurations 



Configuration 


Mean A^p 


Median AMD 


Median RMC 


Median WMF^ 


Median rjo™,gj(Myr) 


CJS 


3.13 


0.010 


29.8 


5.3 X 10-3 


123.5 


CJSECC 


3.5 


0.0047 


40.9 


1.2 X 10-3 


80.4 


EJS 


2.75 


0.0062 


44.2 


3.3 X 10-'' 


147.8 


EEJSi 


3.58 


0.0099 


42.6 


1.1 X 10-* 


109.3 


JSRES 


4 


0.0071 


34.2 


7.5 X 10-3 


160.1 


JSRESECC 


3.75 


0.0041 


28.7 


1.1 X 10-3 


78.6 


Solar System 


4 


0.0018 


89.9 


~ 10-3 


50-150 



"Note that there were 12 EEJS simulations, including 8 with 2000 planetesimals. For the 4 EEJS simulations with 1000 planetesimals 
(EEJS 1-4), the median N, AMD, RMC, WMF^ and Tjorm,® were 3, 0.015, 50.9, 8.1 x 10"^, and 109.3 Myr, respectively. 



3). As expected, the majority of dry planets come from 
the EJS and EEJS samples. 

We define a successful outcome in Fig. [T2] to have 
T/orm.e = 50-150 Myr WMF(s > 5 x 10-"*. The wa- 
ter constraint requires two oceans of water to have been 
accreted by the planet, because Earth's minimum bulk 
water content is two oceans, one on the surface and one 
in the mantle (Lecuyer et al. 1998; 1 ocean — 1.5 x lO^^g 
is the amount of water on Earth's surface). Of our 40 
simulations, 14 were successful in WMF(^ — Tform,® 
space, but there was no overlap with the 4 successful 
cases from RMC - AMD space. The CJS and CJSECC 
cases were the most successful in this respect: 10/12 CJS 
and CJSECC simulations satisfied both the WMF^ and 
Tform,® constraints. 

6. DISCUSSION AND CONCLUSIONS 

In this section we first discuss the degree to which each 
giant planet configuration was able to reproduce our ob- 
served constraints (§6.1). Next, we discuss the context of 
each giant planet configuration in terms of the Solar Sys- 
tem as a whole (§6.2). We then point out the limitations 
of our simulations, and plans for future work (§6.3). 

6.1. Success of giant planet configurations in satisfying 
our constraints 

Let us quantitatively evaluate how the simulations 
fared at reproducing our five constraints using relatively 
generous values: 1) MMars < 0.3 M®, 2) AMD < 0.0036 
(twice the MVEM value), 3) 50 Myr < Tform,® < 150 
Myr, 4) Less than 0.05 M^ in embryos is stranded in the 
asteroid belt, and 5) WMF^ > 5 x 10-*. No single simu- 
lation reproduced all five constraints. Three simulations 
reproduced four constraints. The simulation CJS15-3 re- 
produced four constraints, but formed a 0.98 M^ Mars 
analog at 1.37 AU and a 0.58 M® planet at 1.91 AU. The 
simulation EJSl-3 formed a small Mars, a wet Earth on 
the correct timescale, and stranded one, 0.048 M^ em- 
bryo in the asteroid belt, but the system's AMD is far 
too large (0.012). The simulation EEJS-9 formed a small 
Mars, low-ylAf_D planets and a wet Earth with one 0.04 
M0 embryo in the main belt. Many simulations repro- 
duced three constraints, but rarely while forming a small 
Mars. As expected from previous work. Mars' small size 
was the most difficult constraint to reproduce, and no 
simulations outside of the EJS and EEJS configurations 
form a single Mars analog less massive than 0.5 M0. 

Table 5 crudely summarizes the outcomes of our simu- 
lations in terms of the likelihood of a system with a given 



giant planet configuration's ability to quantitatively re- 
produce our observed constraints using the values above. 
We completed the table as follows. A configuration is 
said to reproduce a given constraint (and receives a "v^") 
if at least half of the simulations were successful for 
that constraint, using the constraints listed immediately 
above. A configuration is unsuccessful at reproducing 
a constraint (and receives a "x") if no simulations are 
successful. If isolated cases or a small fraction of sim- 
ulations are successful, then the configuration receives 
a maybe ("~"). In one case we bent these rules; 5/12 
EEJS simulations formed a Mars analog smaller than 0.3 
M® (4 cases < 0.2 M®) so we gave EEJS a / for this 
constraint despite a slightly less than 50% success rate. 

The most successful giant planet configuration was 
EEJS ("Extra-Eccentric Jupiter and Saturn"). The 
EEJS simulations reliably satisfied three of our con- 
straints with two maybes (see Table 5), and the ensem- 
ble of EEJS simulations satisfied all five constraints, al- 
though no single simulation did so. The EEJS cases re- 
liably formed reasonable Mars analogs in terms of Mars' 
mass, orbit, and formation timescale. Planetary eccen- 
tricities were too large in most cases, but increasing the 
number of planetesimals (runs EEJS 5-12) decreased the 
AMD to close to the MVEM value, and even higher res- 
olution simulations would presumably continue decrease 
the AMD to the MVEM value. The Earth analog formed 
on the correct, 50-150 Myr timescale in most EEJS sim- 
ulations but was too dry in all but three cases. Almost 
half (5/12) of the EEJS simulations finished with no em- 
bryos in the asteroid belt, and for the cases with trapped 
asteroidal embryos they were typically low-mass. How- 
ever, we note that for the EJS and EEJS configurations 
the survival of any embryos in the asteroid belt consti- 
tutes a failure because Jupiter and Saturn's orbits do 
not allow for any migration, which would be necessary 
to clear remaining embryos from the belt - this issue is 
discussed further in §6.2. 

The main reason for the success of the EEJS simula- 
tions was the strength of the secular resonance located 
at ~ 2 AU, which created an "edge" to the inner disk, ef- 
fectively separating it from the asteroid region, removing 
material that approached 2 AU and thereby helping to 
form a small Mars. The high eccentricities of Jupiter and 
Saturn {ejup,sat — 0.07 — 0.10) were responsible for the 
strength of the resonance; the EJS simulations had vq in 
the same place but it was too weak to clear out enough 
material to form a small Mars. However, having the vq 
at this location also makes it difficult for water-bearing 
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TABLE 5 

Summary of the success of Jupiter-Saturn configurations 

FOR reproducing INNER SOLAR SYSTEM CONSTRAINTS"'- 



Config. 


AMD 


^Mars 


^form 


Ast. Belt 




CJS 




X 






/ 


CJSECC 




X 




X 




EJS 












EEJS 


•/ 


^ 








JSRES 


•/ 


X 




X 




JSRESECC 


•/ 


X 




X 





"For each configuration of Jupiter and Saturn, a check ("v^") 
represents success in reproducing a given constraint in at least half 
the simulations, a cross ("x") represents a failure to reproduce the 
constraint in any simulations, and a twiddle sign C'-^") represents 
a "maybe", meaning success in reproducing the constraints in a 
smaller fraction of cases. 

asteroidal material to enter the inner Solar System and 
be accreted by the terrestrial planets. 

In several of the giant planet configurations that we 
considered - CJS, CJSECC, JSRES, and JSRESECC - 
Jupiter and Saturn had lower eccentricities than their 
current values. For all of these simulations, the Earth 
generally formed on the correct timescale and the terres- 
trial planets were low-eccentricity and water-rich. How- 
ever, not a single simulation from these four cases was 
able to reproduce Mars' small size, and most simulations 
stranded one or more large embryos in the asteroid belt. 
Therefore, this work suggests that low-eccentricity con- 
figurations of Jupiter and Saturn cannot explain the ter- 
restrial planets, in particular Mars' small size, in the con- 
text of our simulations. 

It is interesting that none of our simulations was able 
to reproduce the large radial mass concentration seen in 
the Solar System's terrestrial planets {RMC — 89.9 for 
MVEM vs. 30-50 for most simulations; see Fig. [TT] and 
Eqn 1). This concentration comes from the large masses 
and proximity of Venus and Earth, and the small masses 
of Mercury and Mars. Given the difficulty in producing 
Mars analogs, it is not surprising that simulations with 
low-eccentricity giant planets yield small RMC values. 
However, the EEJS and EJS simulations also yielded 
RMC values far smaller than MVEM, although larger 
than for the other giant planet configurations. The ori- 
gin of this discrepancy is not clear. It could be related 
to the structure of the planetesimal disk; observations 
suggest that inner, dust-free cavities exist in many disks 
around young stars with varying radii, from < 0.1 to ^ 
1 AU. (e.g., Eisner et al. 2005, Millan-Gabet et al. 2007). 
If the Solar Nebula had a large inner cavity then the in- 
ner boundary for the planetesimal disk could have been 
at roughly Venus' orbital distance such that the radial 
compression of MVEM is a result of accretion in a radi- 
ally compressed planetesimal disk. Alternately, resonant 
sweeping or tidally-induced migration from interactions 
with the residual gas disk or perhaps coUisional debris 
could compress the terrestrial planet system and increase 
the RMC (see Thommes et al.2008). However, a diffi- 
culty with this model is that Earth's formation timescale 
is much longer than the typical gas disk lifetime (see 
§6.3). 

6.2. Putting the giant planet configurations in the 
context of the Solar System 



Our giant planet configurations, described in §3.1, rep- 
resent different assumptions about the early evolution of 
the Solar System. Four of our cases - CJS, CJSECC, 
JSRES, and JSRESECC - are based on the Nice model, 
which requires Jupiter and Saturn to have formed inte- 
rior to their mutual 2:1 mean motion resonance (Tsiga- 
nis et al. 2005; Morbidelh et al. 2007). These four cases 
assume that delayed, planetesimal scattering-driven mi- 
gration spread out the giant planet system, and that the 
2:1 resonance crossing of Jupiter and Saturn triggered 
the late heavy bombardment (Gomes et al. 2005; see also 
Strom et al. 2005). The other two giant planet config- 
urations - EJS and EEJS - have Jupiter and Saturn at 
their current semimajor axes, meaning that planetesimal- 
scattering driven migration is not permitted because this 
inevitably leads to spreading out of the giant planets' 
orbits (Fernandez & Ip 1984). Therefore, the EJS and 
EEJS simulations assume that the Nice model is incor- 
rect and that the late heavy bombardment was caused 
by another mechanism. The only other currently-viable 
such mechanism is the "planet V" theory of Chambers 
(2007) which invokes the formation and delayed insta- 
bility of a fifth, sub-Mars mass terrestrial planet at ~ 2 
AU. Thus, for EEJS and EJS simulations to be consis- 
tent with the Solar System's observed history they must 
form a planet V with the right mass (< MMars) and in 
the right location (~ 2 AU). 

Our results clearly favor the EEJS simulations over 
giant planet configurations that are consistent with the 
Nice model, mainly because of the EEJS simulations' 
ability to form Mars analogs. In addition, the constraints 
that were poorly reproduced by the EEJS simulations 
were in some sense our weakest, because higher resolu- 
tion simulations with more planetesimals tend to lower 
the AMD, and alternate models exist for water deliv- 
ery to Earth (e.g., Ikoma & Genda 2006; Muralidharan 
et al.2008). Moreover, EEJS simulations 5, 6 and 11 
formed reasonable planet V analogs (see Table 2). In one 
of these cases, simulation EEJS-11, a good Mars analog 
formed at 1.62 AU and a 0.09 M® planet V formed at 2.0 
AU. Given the significant inclination of the planet V in 
this case (12°), the instability timescale could be several 
hundred Myr (e.g., Morbidelh et al. 2001). That simu- 
lation therefore may represent the most self-consistent 
simulation in our sample, at least when only considering 
the inner Solar System constraints. 

Given that planets are thought to form on circular or- 
bits (e.g.. Pollack et al. 1996), what would be the source 
of the significant eccentricities of Jupiter and Saturn in 
the EEJS configuration? Goldreich & Sari (2003) pro- 
posed that eccentricities of Jupiter-mass planets could 
be excited by lindblad resonances in gaseous protoplane- 
tary disks, but with eccentricities limited by the width of 
the gap carved out of the disk by the planet. D'Angelo et 
al. (2006) found that planets embedded in eccentric disks 
can have their eccentricities increased to ~ 0.1. How- 
ever, other studies have shown that eccentricity growth 
via planet-disk interactions can only occur for very mas- 
sive, > IQMjup planets (Papaloizou et al. 2001; Kley & 
Dirksen 2006). Thus, the origin of EEJS-type configu- 
rations of Jupiter and Saturn remains controversial but 
certainly within the realm of current thinking. 

Another study has shown that the EEJS configura- 
tion may be advantageous for reproducing the inner So- 
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lar System. Thommes et al. (2008) also took advan- 
tage of strong secular resonances to reproduce the ter- 
restrial planets. In their model, collisions between em- 
bryos were induced by the inward sweeping of the 
and uq, which occurred as the Solar Nebula dissipates 
and changes the gravitational potential (Heppenheimer 
1980; Ward 1981). Their source of damping is tidal in- 
teraction with the gas (Ward 1993; Cresswell et al. 2007). 
They manage to reproduce several aspects of the terres- 
trial planets, including Mars' small size. However, we 
note that Thommes et al. (2008) assumed that the inner 
disk only contained embryos out to 3 AU and no plan- 
etesimals. 

Given these successes, it is tempting to regard the 
EEJS configuration as the true configuration of Jupiter 
and Saturn early in Solar System history. EEJS is indeed 
consistent with the inner Solar System. In the EEJS sce- 
nario, Jupiter and Saturn must have acquired eccentric- 
ities of ~ 0.1 at an early stage, within 1 Myr or so af- 
ter their formation. Scattering of unstable planetesimals 
and embryos over the ~ 10^ years of terrestrial accretion 
decreased their eccentricities to their current values of 
~ 0.05. Indeed, the final time-averaged values of ejup 
and esat in our simulations are 0.03-0.06 and 0.06-0.10, 
respectively, very close to their current orbits. In this 
model. Earth's water was delivered in part from hydrated 
asteroidal material, but mainly from adsorption of small 
silicate grains (Muralidharan et al. 2008), cometary im- 
pacts (Owen & Bar-nun 1995) or oxidation of a H-rich 
primitive atmosphere (Ikoma & Genda 2006). The late 
heavy bombardment of the terrestrial planets can then 
be explained by the delayed destabilization of planet V 
(Ghambers 2007), analogs of which did indeed form in 
several EEJS simulations. 

The outer Solar System provides a strong argument 
against the EEJS configuration. The resonant structure 
of the Kuiper Belt requires outward migration of Nep- 
tune (Malhotra 1993, 1995; Levison & MorbideUi 2003). 
It is thought that Neptune migrated outward because of 
the back-reaction from the scattering of a many Earth 
masses worth of remnant planetesimals. Given that 
Neptune cannot easily eject these planetesimals (e.g., 
Duncan et al. 1987), the scattering of these small bod- 
ies also causes Uramis and Saturn to migrate outward, 
and Jupiter, which does eject most of the small bod- 
ies, to migrate inward (Fernandez & Ip 1984; Hahn & 
Malhotra 1999). This planetesimal scattering is thought 
to have started during planet formation and lasted for 
10^ — 10*^ years, with more rapid migration correspond- 
ing to a shorter migration, a more massive planetesimal 
disk or closer proximity between Neptune and the disk's 
inner edge (Gomes et al. 2004, 2005). Excess depletion 
just exterior to strong resonances in the asteroid belt 
provides empirical corroboration of giant planet migra- 
tion (Minton & Malhotra 2009), although the migration 
timescale cannot be constrained. Thus, the current or- 
bits of Jupiter and Saturn are not thought to be their 
orbits at the time of formation, and it is their orbits at 
early times that affected terrestrial planet formation. In 
particular, Jupiter and Saturn must have formed in a 
more compact configuration, because scattering-induced 
migration always causes their orbits to diverge. 

Is it possible to reconcile the EJS and EEJS configu- 
rations with the outward migration of Neptune that is 



required by current Kuiper Belt models? If a small Mars 
formed because of the resonance, then the t'g miist 
have been in roughly its current location during Mars' 
formation. Hf/W isotopes suggest that Mars formed 
much faster than Earth, and did not undergo any giant 
impacts after 1-10 Myr (Nimmo & Kleine 2007). In ad- 
dition, the location of the is highly sensitive to Jupiter 
and Saturn's semimajor axes such that a more compact 
configuration places the i/g in the asteroid belt, farther 
from Mars and less likely to influence its growth (see Fig. 
2 of Minton & Malhotra 2009). Migration of Jupiter and 
Saturn must therefore have been completely flnished in a 
few Myr at the latest to form a small Mars. Would such 
an early migration fit our current understanding? At 
such early times there likely existed a population of em- 
bryos in the asteroid belt (Wetherill 1992). If Jupiter 
and Saturn's migration occurred that early, then em- 
bryos would probably have smeared out the observed de- 
pletion of asteroids exterior to asteroid belt resonances 
(Minton & Malhotra 2009). In addition, this very rapid 
migration must also have taken place in the presence of 
some amoimt of residual disk gas, which would certainly 
have affected the scattering dynamics (Gapobianco et 
al. 2008), and would likely have been far too efficient in 
trapping Kuiper Belt objects in resonances with Neptune 
(e.g., Mandell et al.2007). This suggests that the EEJS 
configuration is at best marginally consistent with the 
outward migration of Neptune needed by models of the 
Kuiper Belt. Neptune's migration would have to have oc- 
curred early and been very rapid. In addition, to remain 
consistent with the late heavy bombardment, a fifth ter- 
restrial planet must have formed at 2 AU. (Chambers 
2007). 

An alternate hypothesis for the evolution of the outer 
Solar System suggests that there initially existed several 
additional ice giant planets (Chiang et al.2007). Could 
such a scenario be consistent with the EEJS configura- 
tion? In this model, instabilities among the ice giants 
would have led to the ejection of excess ice giants and dy- 
namical friction could re-circularize the orbits of Uranus 
and Neptune (Ford & Chiang 2007). However, in order 
to avoid disrupting the EEJS configuration by injecting 
planetesimals into the Jupiter-Saturn region, those plan- 
etesimals would need to have been dynamically ejected 
by the ice giants themselves, which is unlikely given their 
relatively small masses. In addition, a detailed simula- 
tion of this model suggests that a population of many 
ice giants would not eject each other but would sim- 
ply spread out by interactions with the planetesimal disk 
(Levison & MorbideUi 2007). Thus, we are left with the 
same issue, that the giant planets would have had to clear 
out the outer Solar System in less than Mars' formation 
timescale of a few Myr. This scenario therefore does not 
appear to be consistent with the EEJS configuration. 

The Oort cloud provides an additional argument 
against the EEJS configuration. If the giant planets 
reached their current orbits very quickly (as required by 
EEJS), then the Oort cloud had to form very quickly 
as well. To populate the inner and outer classical Oort 
cloud, the current galactic environment (or one moder- 
ately denser) is required (Brasser et al. 2006; Kaib & 
Quinn 2008). However, to capture Sedna (Brown et 
al. 2004), a significantly denser environment is needed. 
These two constraints are not at odds if the giant plan- 
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ets evolved on long timescales, because the galactic en- 
vironment can change on a 10-100 Myr timescale (e.g., 
Lamers et al. 2005). However, if the giant planets finished 
migrating, and therefore finished clearing the planetesi- 
mal disk, within a few Myr, then the Solar System should 
have either an Oort cloud or Sedna, but not both. 

Thus, we are left with a problem: the giant planet con- 
figuration which best reproduces the terrestrial planets 
(EEJS) is at best marginally consistent with the current 
view of the outer Solar System's evolution, and likely in- 
consistent. On the other hand, configurations that are 
based on current outer planet evolution models, in par- 
ticular the Nice model, cannot form a small Mars, an 
inescapable inner Solar System constraint. Given that 
we do not have a choice for the configuration of Jupiter 
and Saturn at early times that satisfies all of our con- 
straints and is also consistent with the evolution of the 
outer Solar System, we cannot reject the Nice Model. 

6.3. Simulation limitations and future work 

It is important to note that our simulations are missing 
several physical effects that could be important. These 
effects include collisional fragmentation (Alexander & 
Agnor 1998), dynamical effects of collisional debris (Lev- 
ison et al. 2005), tidal damping and resonant sweeping 
from a residual amount of nebular gas (Kominami & Ida 
2002; Nagasawa et al. 2005; Thommes et al. 2008). 

We do not think that including collisional fragmenta- 
tion would affect our results for two reasons: 1) Alexan- 
der & Agnor (1998) saw little difference in their simula- 
tions when they included a simple fragmentation model, 
and 2) the velocities of embryo-embryo collisions were 
slow enough that the majority were accretionary rather 
than erosive (Agnor & Asphaug 2004). 

We also assumed that all embryo-planetesimal impacts 
were accretionary, although these do tend to occur at 
larger velocities than embryo-embryo impacts because 
of viscous stirring (see Fig. [5]). Melosh (1984) showed 
that it is difficult for km-sized impactors to accelerate 
material to speeds greater than one third of the impact 
speed. Based on this, we can assume that planetesimal 
impacts at speeds of less than 3 times the escape speed 
Vesc are 100% accretionary; 83% of the collisions in our 
simulations meet this criterion. Furthermore, only for 
impacts with v/vesc > 7 — 10 do small-body collisions ac- 
tually become erosive (O'Keefe & Ahrens 1977; Svetsov 
2007). Only 2 out of 3400 collisions in our simulations 
had V I Vesc > 7. Thus, our assumption of perfect accre- 
tion for planetesimal-embryo impacts probably does not 
affect our results if planetesimals were indeed km-sized. 

The dynamical effects of collisional debris are certainly 
important, although their effects have barely been ex- 
plored because of numerical limitations. A continuous 
source of small bodies after large impacts would likely 
reduce the eccentricities of the terrestrial planets and 
make it easy to reproduce the low AMD of MVEM (Lev- 
ison et al. 2005). In addition, planetesimal-planetesimal 
collisions create collisional cascades that can grind plan- 
etesimals to dust, which can then be removed from the 
system (e.g., Kenyon & Bromley 2006). However, inter- 
actions between the gas disk and small fragments at early 
times may accelerate embryo growth by stopping colli- 
sional cascades and allowing embryos to accrete small 
fragments (Kenyon & Bromley 2009). 



Tidal damping from small amounts of remnant disk gas 
can help reduce the eccentricities of the terrestrial plan- 
ets, as shown by Kominami & Ida (2002, 2004) and Ag- 
nor & Ward (2002). Gaseous disks around other stars are 
thought to disperse in 5 Myr or less (Haisch et al. 2001; 
Briceno et al. 2001; Pascucci et al. 2006). Such lifetimes 
correspond to a decrease in gas density to a level of about 
10% of the minimum-mass nebular density, while tidal 
damping can operate down to a level of 10~^ to 10""*. 
Nonetheless, for this to be important the gas density 
must be at least 10~^ for the timescale of Earth's ac- 
cretion, 50-150 Myr (Touboul et al. 2007). For a steady 
accretion model, the gas density E decreases relatively 
slowly with time t, as t~^^^ (Lynden-Bell & Pringle 
1974). Thus, if the Solar Nebula evolved smoothly with a 
characteristic timescale of 1 Myr, then its density would 
have decreased to 10% after 5 Myr, but would still be 
10-3 after 100 Myr. In such a scenario, tidal damping 
would indeed be an important effect on terrestrial ac- 
cretion. However, the final phases of disk dissipation are 
thought to occur after a few Myr (Haisch et al. 2001) on a 
~ 10^ year timescale (Simon & Prato 1995; Wolk & Wal- 
ter 1996) and relatively violently, via photo-evaporation 
(HoUenbach et al. 1994; Johnstone et al. 1998) or poten- 
tially the MRI instability (Chiang & Murray-Clay 2008). 
Future, more sensitive observations that probe smaller 
gas densities in disks around young stars will shed light 
on this issue, but our interpretation of the current state 
of knowledge is that it is unlikely that tidal damping 
would dramatically change our results. Indeed, we sus- 
pect that the most important source of damping during 
accretion is likely to be small bodies, i.e., planetesimals 
and collisional debris. 

Secular resonance sweeping during the dispersal of the 
Solar Nebula could play an important role in terres- 
trial planet formation if Jupiter and Saturn's orbits have 
eccentricities of at least ~ 0.05 (Nagasawa et al. 2005; 
Thommes et al. 2008). However, if the giant planets' or- 
bits are less eccentric then secular resonances are too 
weak to have much effect on the outcome (O'Brien et 
al. 2007). We note that the model of Thommes et 
al. (2008) invokes this resonant sweeping to induce col- 
lisions as well as to shepherd embryos from the aster- 
oid belt inward. However, Thommes et al. assume that 
the entire inner disk is composed of embryos and do not 
include any planetesimals in their simulations. Given 
the strength of aerodynamic gas drag on km-sized bod- 
ies (Adachi et al. 1976), planetesimals should also have 
been shepherded inward by the resonance (e.g., Raymond 
et al. 2006b; Mandell et al. 2007). Indeed, the mass dis- 
tribution in the area swept out by the secular resonances 
may have been drastically altered if the resonances were 
strong enough. We are currently exploring the conse- 
quences of this idea. 

We did not test an exhaustive number of configurations 
of Jupiter and Saturn. Indeed, given the relative suc- 
cess of the EEJS configurations, it would be interesting 
to explore the resonant behavior for other allowed con- 
figurations such as, for example, Jupiter and Saturn on 
their CJS orbits but with eccentricities of 0.1. We plan 
to search the parameter space of giant planet configu- 
rations in future work, focusing on systems with strong 
mean motion and secular resonances in the inner Solar 
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System. 

The majority of this paper has dealt with the effects 
of the giant planet configuration on the accretion of the 
inner Solar System. However, we note that the disk's 
surface density distribution could also play a strong role 
in the outcome (Raymond et al. 2005). Indeed, we did see 
changes in the 8 simulations we performed with a flatter, 
surface density profile. We note that there exists an 
alternate model that attributes Mars' small mass with 
the Solar Nebula's density structure (Jin et al. 2008). In 
Jin et al. 's model the disk is ionized only in certain re- 
gions, causing radial variations in the strength of the 
magneto-rotational instability and therefore in the disk 
viscosity (Balbus & Hawley 1991). At the interface be- 
tween an outer low- viscosity and an inner high-viscosity 
regime a local dip in the surface density can be created. 
Jin et al. (2008) 's model has this dip at about Mars' or- 
bital distance. This dip is quite deep but very narrow, 
although it could have swept over a region of radial with 
~1 AU in the lifetime of the disk. When considering 
accretion in such a disk, it is important to note that 
the typical planetary feeding zone in our sample has a 
width of 2-3 AU, much wider than the widest possible 
gap created in the disk of Jin et al. . Indeed, preliminary 
simulations using Jin et al. 's disk form Mars analogs that 
are ~ 5MMars- 

Our simulations are among the highest-resolution to 
date, but we are still resolution-limited. To accurately 
model Mars' growth we require that the embryo mass 
at Mars' orbital distance is smaller than Mars itself, in 
our case by a factor of 3-6 (Fig. [J). This mass sets the 
inter-embryo spacing and we are subsequently limited by 
the number of particles we can integrate in a reasonable 
time (of several months). Given these restrictions, our 
simulations may not always adequately model dynamical 
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friction, because the planetesimal-to-embryo mass ratio 
is relatively small (~ 10 at Mars' orbital distance; smaller 
closer in and larger farther out). Indeed, the EEJS sim- 
ulations with 2000 planetesimals yielded smaller eccen- 
tricities than those with 1000 planetesimals because we 
could resolve the damping at late times. Our approach 
contrasts with that of Chambers (2001) and O'Brien et 
al. (2006), who used fixed-mass embryos and therefore 
had fixed planetesimal-to-embryo mass ratios of 10 and 
40, respectively. The advantage of their approach is that 
dynamical friction is more consistently modeled, but the 
disadvantage is that the number of embryos is small and 
the embryo mass is large, about one Mars mass. Despite 
the limited resolution and the difficulty with dynamical 
friction, all of our constraints were met in some simula- 
tions, including cases with AAID values lower than the 
MVEM value. We therefore think that our approach is 
valid, and we anticipate that faster computers will cer- 
tainly improve the outlook for understanding the origin 
of the inner Solar System in the coming years. 
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